A robust upper limit on N e g from BBN, circa 2011 
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We derive here a robust bound on the effective number of neutrinos from constraints on primordial 
nucleosynthesis yields of deuterium and helium. In particular, our results are based on very weak 
assumptions on the astrophysical determination of the helium abundance, namely that the minimum 
effect of stellar processing is to keep constant (rather than increase, as expected) the helium content 
of a low-metallicity gas. Using the results of a recent analysis of extragalactic HII regions as upper 
limit, we find that AN c g < 1 at 95% C.L., quite independently of measurements on the baryon 
density from cosmic microwave background anisotropy data and of the neutron lifetime input. In 
our approach, we also find that primordial nucleosynthesis alone has no significant preference for an 
effective number of neutrinos larger than the standard value. The ~ 2<r hint sometimes reported 
in the literature is thus driven by CMB data alone and/or is the result of a questionable regression 
protocol to infer a measurement of primordial helium abundance. 
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I. INTRODUCTION 



Historically, the helium abundance has played an im- 
portant role for establishing the "hot big bang" cosmo- 
logical model, see for example [lj . Primordial nucleosyn- 
thesis, most often referred to as big bang nucleosynthesis 
(BBN), helped building confidence in the overall credi- 
bility of cosmology as a science. Nowadays, in the era of 
"precision cosmology" , it is clearly no more question to 
prove that the bulk of 4 He is primordial, also in view of 
the indirect (but clean) detection of a non- vanishing 4 He 
mass fraction Y p from CMB data, see 0, Q . 

The basic pillars of the hot big bang model have found 
many confirmations and have been subject to several 
important cross-checks, leading to the so-called concor- 
dance (or standard) model of cosmology. In the last two 
decades, BBN has thus mostly turned into a probe of 
physics/cosmology beyond the standard model (see [3,0 
for recent reviews). Moving forward in this direction 
has been however hampered by the systematics in the 
determination of light nuclei abundances, since one has 
no access to truly primordial environments, i.e. prior to 
the first generation of stars and nucleosynthetic events: 
any inference thus relies to some extent on astrophys- 
ical models. In particular for the case of 4 He, despite 
the more and more careful analyses undertaken over the 
past decade, very different conclusions can be sometimes 
found on the inferred Y p , see e.g. jg-|l|. Without en- 
tering the issues related to the analysis of spectroscopic 
data, from a particle physics perspective one may wonder 
what is the most profitable way to use the results of these 
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analyses for constraining new physics, while maintaining 
some robustness and independence from the regression 
protocol. 

In fact, we believe that for a largest fraction of the par- 
ticle astrophysics community it is more important to ob- 
tain reliable constraints than inferring primordial abun- 
dances. In this spirit, here we propose a simple and ro- 
bust approach to the use of 4 He data for constraining 
the effective number of neutrinos, N c s, by far the most 
widely used BBN-related quantity used to parameterize 
new physics. This Letter is structured as follows: in 
Sec. |TT] we outline our minimal assumptions in using the 
data; in Sec. IIIII we present our results; in Sec. IIVI we 
discuss our findings and assumptions and conclude. 



II. PROCEDURE 

For 4 He, we use here the abundances inferred in nine 
metal-poor, extragalactic HII regions in Q (see also 0). 
Differently from the usual practice, we do not perform a 
regression to zero metallicity, since our aim is to derive 
an upper bound to 4 He. Hence, we fit the yields to a 
constant abundance Yq, obtaining 

(Y ) ± era = 0.2581 ± 0.0025 (68% C.L.) . (1) 

which also implies 

Y < 0.2631 at 95% C.L. (2) 

Assuming that no information is available on the lower- 
limit of 4 He, one can parametrize the likelihood function 
for Y p as a flat+semi-gaussian shape 



'Electronic address: mangano@na.infn.it 



£(Y P ) oc Q((Y )-Y p )+Q(Y p -<Xo))exp 



^Electronic address: serpico@lapp.in2p3.fr 



(Y P -(Yo)) 2 
2al 



(3) 



2 



We argue here that this is conservative and robust. In 
fact, Eq. @ is a strict upper limit to the primordial value 
Y p under the sole assumption that dY/dZ > 0, i.e. that 
the average 1 effect of stellar processing is to increase 
the helium content of a low-metallicity gas. By now, 
there is some empirical evidence of this trend (positive 
derivative) even in the few high-quality objects analyzed 
in 0, Hj]. However, assuming dY/dZ =constant is a much 
stronger assumption, usually justified empirically since 
the seventies as the simplest possible fit [|| (see also [Io|). 
but subject to the obvious risk of extrapolation errors. 

Even if one could model the metallicity evolution of 
the observed systems reliably, it is unclear to what ex- 
tent the pre-galactic value of Y should coincide with its 
primordial value Y p . Actually, it has been proven since 
longtime that the bulk of 4 Hc must be primordial (see 
e.g. [ill ] for a short historical overview), and by now there 
is a positive (albeit indirect) detection of a non- vanishing 
Y p from CMB data, see 0, H[- However, it is a differ- 
ent issue to prove that no extra production of 4 He has 
taken place since early times at a level comparable to the 
current error on its determination, i.e. at the percent 
level. In fact, "What is the helium enrichment by the 
first stellar generations?" was a question already posed 
in the seminal paper In the last decade, several 

models have been proposed where such a production does 
take place at an appreciable level (AY ~ 10~ 3 — 10 -2 ), 
see for example (l3l flij . This is perhaps not surpris- 
ing, given that the yet undiscovered generation of stars 
known as PopIII forming in the pristine metal-free gas 
should generate some chemical enrichment. It is worth 
remembering that currently observed "metal-poor" sam- 
ples have metal contents many orders of magnitude above 
the expectations from BBN [15]. When using the late 
universe determinations of Y as measurements of its pri- 
mordial value, a systematics error is committed, difficult 
to quantify precisely since it depends on astrophysical 
modeling of pre-galactic times. Again, this extra prob- 
lem is avoided in our conservative procedure. Notice that 
our upper bound of Eq. ([2]) is close but slightly more 
stringent than what is found at 2 a in H, Y p < 0.2639, 
since only seven of the nine determinations were used 
there. The small difference shows however that a mi- 
nor change in the choice of objects would not change the 
following results appreciably, which is another hint in fa- 
vor of their robustness. It is more involved to compare 
with the results of the group of Izotov et al., which does 
not present fits assuming dY/dZ = 0. For illustration, 
the measurement quoted in [6j] , based on a linear regres- 
sion method whose robustness we question here, yielded 
Y p = 0.2565 ± 0.0010 (la stat.) ± 0.0050 (syst.). On the 
other hand, we note that fitting to a constant the 10 
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determinations reported in Table 4 of [l6j], as suggested 
here, one obtains Y = 0.2555 ± 0.0016, i.e. Y < 0.2587 
at 95% C.L. When rescaling the above value upwards by 
2% due to some improved atomic corrections as argued 
in Q, one finds a value remarkably consistent with the 
above determinations of the upper limit. We take this 
exercise as an indication that our results are not sensi- 
tive to the different analysis codes and protocols, at least 
when the same and most updated atomic data are used. 

For 2 H, following the analysis in [4], as best estimate 
of the primordial deuterium yield one may use 

( 2 H/H) P = 2.87±°; 22 at 68% C.L. (4) 

which is based on a conservative analysis of seven (rela- 
tively) reliable absorption spectra of clouds at high red- 
shifts, on the light of background quasars. Detailed stud- 
ies [10| suggest that the depletion due to stellar activity 
in these early systems is negligible and thus the above 
range can be considered a faithful estimate of the primor- 
dial value. Yet, we shall also show one example based on 
the minimal assumption 

( 2 H/H) P > 2.45 x 10~ 5 at 95% C.L. , (5) 

which is agnostic on an eventual depletion of 2 H in an 
early chemical evolutionary phase. In this case, a semi- 
gaussian+fiat likelihood function similar to what done 
for Helium in Eq. ([3]) is constructed. Even in this case, 
very similar results follow on the upper limit on N c g. 

Concerning the BBN predictions, we used nuclide 
yields based on the code described in [l7|, also includ- 
ing the uncertainties due to nuclear reactions extensively 
described in [l8j |. 

Finally, one may (or may not) impose the CMB mea- 
surement of the baryon fraction [2] 

u b = n b h 2 = 0.02250 ± 0.00056 at 68% C.L. (6) 

We checked that significantly relaxing the above range 
(say, by one order of magnitude) does not affect the upper 
limit on N e g appreciably. Note that for consistency we 
consider the value for uib inferred in a ACDM+iV c ff model 
(i.e. where N c g is allowed to vary), although using the 
best-fit in a ACDM model would make no quantitative 
difference. 

For illustration, in one case we shall also consider the 
impact of the Y p measurement from CMB data (WMAP 
7 yrs plus small scale results of the Atacama Cosmol- 
ogy Telescope Q ) , assuming the additional measurement 
(gaussian error) 

Y p = 0.313 ± 0.044 at 68% C.L. (7) 

III. RESULTS 

With the ingredients described above we constructed 
two-dimensional likelihood functions in the Wb — N e Q 
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TABLE I: Constraints on N e g corresponding to different 
datasets used: i) first row: Eq. ((6|, Eqs. (|1I3[1 . Eq. (j4]); ii) sec- 
ond row: Eq. ©, Eqs. PT3)l, Eq. ©; iii) third row: Eqs. (pT3)) 
and Eq. Q; fourth row: as the first one, with the additional 
CMB measurement of Y p of Eq. ©. The last column shows 
the likelihood that N b b is smaller than the standard value 
3.046 [JJ]. 

plane, then marginalized over the parameter u>b, which 
is not of interest here. The results of our analysis are 
thus encoded in the 1-dimensional likelihoood functions 
£(iVeff), whose integrals are normalized to 1. These func- 
tions are shown in Fig. [TJ and relevant numerical quan- 
tities are summarized in Table [TJ We define N^ ln and 
A^ ax such that 

/ C(x)dx = 0.95 , / C(x)dx = 0.95, (8) 
ijv™ f in Jo 

and the parameter L(N eS < N^) in Table Has 

L(N ee <N™)= / C(x)dx. (9) 
Jo 

When remembering that the standard model expec- 
tation for N e g is about 3.046 we see that in all 
cases we get a bound A7V cff < 1. The reason why it is 
slightly more stringent when using deuterium as a "mea- 
surement" of cj;, instead of CMB (third line) is that it 
favors a slightly smaller value for the baryon fraction. 
In correspondence of this smaller value, the deuterium 
yield is a bit larger. Since deuterium grows with N c g, in 
this case the deuterium hits the upper bound for a lower 
value of N c g, increasing its constraining power. As it is 
clear from the second row of Table [TJ allowing for primor- 
dial deuterium depletion and limiting oneself to consider 
the lowest limit of its measured value as a lower limit, 
the bound does not change much, since the constraining 
power derives from the upper limit on He. In Fig. [JJ 
this reflects on the quite hard cut in the likelihood func- 
tions at large N c g. Also, adding the CMB measurement 
of Y p of Eq. ([7]) does not change much the situation with 
respect to the first case: the slight shift towards higher 
values of N e g reported in the fourth row is simply due to 
the fact that the current best value of Y p from CMB is 
above the BBN prediction, albeit not significantly (less 
than 1.5 a). This also proves indirectly that if we had 
imposed a loose lower-bound on Y p (say, Y p > 0.225) in- 
stead of the flat likelihood of Eq. ([3]) at low- Yp, the result 
would hardly change. 

On the other hand, comparing the first and last two 
lines in the table shows that an independent constraint 



1.0 





1 1 1 1 1 1 1 


i i 




w L +"H+ 4 He 

b : 






2 4 : 

K> b + Hj + He : 


f \ 




2 H+ 4 He 


I \ 




■ ■ ■ Cft+Y + H+ He ■ 

b : 












' /." 

y' /.' 


\ \ 
\\ 


■1-flJ Lid I.I 1 WI-fPT* 1 


\ 

. ■ , , , l 



0.0 1.0 2.0 3.0 4.0 5.0 

N eff 



FIG. 1: Marginalized 1-D likelihood functions C versus N e g 
using the different combinations of data as in Table [I] Solid 
(red) and dashed (purple) curves are obtained using CMB 
measurement of cot, with the dotted (black) one also adds 
CMB information on Y p . In all cases the quite sharp cut-off 
at N e s ~ 4 is due to 4 He abundance upper limit. 

on Ldf, and possibly even a relatively weak lower limit on 
Y p are quite useful in setting a stringent lower limit on 
N e g (second column of Table fl]). In particular, the effect 
of the constraint on uji, is explained as follows: since the 
dependence of He on w;, is very weak, and 2 H suffers 
of a partial degeneracy between N e g and w&, relatively 
low values of N e g can be compensated with relatively 
high values of Ub- Hence, imposing an upper limit on ujb 
yields to a more stringent lower limit on N e g. Of course, 
this exercise has only illustrative purpose: the physics 
behind the CMB measurement on uii, is well understood, 
and any cosmologically meaningful lower limit on N e g is 
significantly larger than the value reported at the third 
row in Table [TJ 

Finally, it is interesting to consider the last column in 
Table [ij illustrating the likelihood that the inferred N e g 
value is lower or equal than its standard model expec- 
tation: we see that BBN alone has no clear preference 
for a larger-than-standard N e g (compared to a lower- 
than-standard one) when the observed abundances are 
interpreted conservatively. The blue, dot-dashed curve 
in Fig. [TJ also shows graphically the same effect. Even 
when combined with CMB data, BBN does not favor 
significantly larger-than-standard values for N e g. 

IV. DISCUSSION AND CONCLUSIONS 

In this letter, we have discussed a new and more con- 
servative approach to derive BBN constraints on N e g, 
motivated by growing concerns on the reliability of astro- 
physical determinations of primordial 4 He. We showed 
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that even with present data, an analysis which is at the 
same time conservative and informative is possible, in 
particular when using CMB and/or deuterium measure- 
ments as constraints on the baryon abundance. In all 
the cases we considered, we found that AN e g < 1 at 95% 
C.L., and possibly the bounds are slightly more stringent, 
see Table |U It is perhaps useful to note that, following 
hints from neutrino laboratory experiments, models with 
two sterile neutrinos with relatively large mixings have 
been recentl y ( re) considered in a cosmological context, 
see e.g. (20l l2l|. Although CMB data alone might ac- 
commodate or even marginally favor such amount of ex- 
tra radiation, our results suggest that they are strongly 
disfavored or excluded by BBN, at least if the two sterile 
neutrino distributions are close enough to thermal ones. 
Of course, this is only an example: non-standard values 
of N e g (and in particular low ones) might even indicate 
new physics completely unrelated to neutrinos (more de- 
tails can be found e.g. in [4J). 

Let us briefly discuss the sensitivity of our results to 
some underlying assumptions. Let's assume that our sole 
astrophysical assumption on helium, dY/dZ > 0, should 
be empirically proven wrong. In such a case (at the mo- 
ment of academic interest, given the evidence of the con- 
trary from the data!), lacking a model-independent way 
to correct for the evolutionary effects, the whole logic 
of "primordial helium traceability" would be put under 
discussion: otherwise said, we expect that observational 
progresses will either improve over our upper bound (for 
example by limiting the analysis to very low metallic- 
ity objects and/or increasing the statistics), or one will 
have to give up completely the idea that any quantita- 
tive inference on the primordial value can be done from 
what observed in low-redshift stellar objects. Only indi- 
rect probes as the CMB one could then be used reliably, 
though with quite a large uncertainty. 

On the other hand, albeit only in model-dependent 
frameworks, one has an idea of the possible contribu- 
tions of pre-galact ic g enerations of stars to Y. For exam- 
ple, the model in [131 ] suggest possible production yields 
of PopIII stars up to AY ~ 3.3 x 10~ 3 , with compara- 
ble effects found in From a numerical evaluation 
of the derivative of Y p (oj£ MB , N c g) at the standard value 
we find AN e g ~ 75 AY, hence a contribution like the 
one above would result in an apparent rightward shift of 
AiVcff ~ 0.25 in the curves of Fig. 1. In [lj] the overall 
chemical evolution is found to be responsible for shifts 



up to AY ~ 0.007 (i.e. AN oS > 0.5) while the impact on 
Deuterium abundances in high-redshift objects is negli- 
gible. It appears reasonable to conclude that the cosmo- 
logical evolution Y(z) due to both pre-galactic and galac- 
tic populations of objects might easily explain the slight 
preference for larger-than-standard N e g values from the 
BBN analysis based on current astrophysical samples. 

One may also wonder how sensitive this bound is to 
other particle physics parameters entering the 4 He pre- 
diction. The most important such parameter is proba- 
bly the neutron lifetime, which sets the timescale of the 
neutron/proton density ratio freeze-out. The two recent 
determinations [22|, H3] of the neutron lifetime differ for a 
5.6 and 2.5 cr's respectively, from the PDG recommended 
world average [24[ used in deriving our results above. Yet, 
we find that adopting the new determinations of the neu- 
tron lifetime would only relax the bound by AN c g <0.1. 

As already noted in the past, the synergy of BBN and 
CMB is useful in narrowing some parameter uncertain- 
ties, once again even under very weak/conservative as- 
sumptions on the measured helium mass fraction. Forth- 
coming data from the Planck experiment are expected 
to strongly increase the sensitivity to AiV e ff, while be- 
ing insensitive to astrophysical pollution of primordial 
Y p . Values for this parameter compatible with BBN 
data expectations, as for example those presented in this 
letter, could then be further refined by exploiting the 
theoretical known dependence of the 4 He mass fraction 
Y p (uJb, N c s) in the analysis, thus reducing the impact of 
a flat prior on Y p on cosmological parameter estimation, 
see e.g. [25|,[26|]. On the other hand, evidence for larger 
than standard N e g from CMB would be a real break- 
through, strongly supporting non-minimal scenarios, as 
the presence of large neutrino chemical potentials [271 ] . 
or — in case of very large values of N e g — hinting for extra 
species changing the evolution of the Universe possibly 
after BBN, but before/during CMB formation, one re- 
cent example of which has been provided in [28]. 
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